Distinct functions of POT1 at telomeres.

The mammalian protein POT1 binds to telomeric single-stranded DNA (ssDNA), protecting chromosome ends from being detected as sites of DNA damage. POT1 is composed of an N-terminal ssDNA-binding domain and a C-terminal protein interaction domain. With regard to the latter, POT1 heterodimerizes with the protein TPP1 to foster binding to telomeric ssDNA in vitro and binds the telomeric double-stranded-DNA-binding protein TRF2. We sought to determine which of these functions-ssDNA, TPP1, or TRF2 binding-was required to protect chromosome ends from being detected as DNA damage. Using separation-of-function POT1 mutants deficient in one of these three activities, we found that binding to TRF2 is dispensable for protecting telomeres but fosters robust loading of POT1 onto telomeric chromatin. Furthermore, we found that the telomeric ssDNA-binding activity and binding to TPP1 are required in cis for POT1 to protect telomeres. Mechanistically, binding of POT1 to telomeric ssDNA and association with TPP1 inhibit the localization of RPA, which can function as a DNA damage sensor, to telomeres.

Interaction of Cryptococcus neoformans Rim101 and protein kinase A regulates capsule.

Cryptococcus neoformans is a prevalent human fungal pathogen that must survive within various tissues in order to establish a human infection. We have identified the C. neoformans Rim101 transcription factor, a highly conserved pH-response regulator in many fungal species. The rim101 multiply sign in circle mutant strain displays growth defects similar to other fungal species in the presence of alkaline pH, increased salt concentrations, and iron limitation. However, the rim101 multiply sign in circle strain is also characterized by a striking defect in capsule, an important virulence-associated phenotype. This capsular defect is likely due to alterations in polysaccharide attachment to the cell surface, not in polysaccharide biosynthesis. In contrast to many other C. neoformans capsule-defective strains, the rim101 multiply sign in circle mutant is hypervirulent in animal models of cryptococcosis. Whereas Rim101 activation in other fungal species occurs through the conserved Rim pathway, we demonstrate that C. neoformans Rim101 is also activated by the cAMP/PKA pathway. We report here that C. neoformans uses PKA and the Rim pathway to regulate the localization, activation, and processing of the Rim101 transcription factor. We also demonstrate specific host-relevant activating conditions for Rim101 cleavage, showing that C. neoformans has co-opted conserved signaling pathways to respond to the specific niche within the infected host. These results establish a novel mechanism for Rim101 activation and the integration of two conserved signaling cascades in response to host environmental conditions.

Mass spectrometry-based thermal shift assay for protein-ligand binding analysis.

Described here is a mass spectrometry-based screening assay for the detection of protein-ligand binding interactions in multicomponent protein mixtures. The assay utilizes an oxidation labeling protocol that involves using hydrogen peroxide to selectively oxidize methionine residues in proteins in order to probe the solvent accessibility of these residues as a function of temperature. The extent to which methionine residues in a protein are oxidized after specified reaction times at a range of temperatures is determined in a MALDI analysis of the intact proteins and/or an LC-MS analysis of tryptic peptide fragments generated after the oxidation reaction is quenched. Ultimately, the mass spectral data is used to construct thermal denaturation curves for the detected proteins. In this proof-of-principle work, the protocol is applied to a four-protein model mixture comprised of ubiquitin, ribonuclease A (RNaseA), cyclophilin A (CypA), and bovine carbonic anhydrase II (BCAII). The new protocol's ability to detect protein-ligand binding interactions by comparing thermal denaturation data obtained in the absence and in the presence of ligand is demonstrated using cyclosporin A (CsA) as a test ligand. The known binding interaction between CsA and CypA was detected using both the MALDI- and LC-MS-based readouts described here.

Thermodynamic analysis of ligand-induced changes in protein thermal unfolding applied to high-throughput determination of ligand affinities with extrinsic fluorescent dyes.

The quantification of protein-ligand interactions is essential for systems biology, drug discovery, and bioengineering. Ligand-induced changes in protein thermal stability provide a general, quantifiable signature of binding and may be monitored with dyes such as Sypro Orange (SO), which increase their fluorescence emission intensities upon interaction with the unfolded protein. This method is an experimentally straightforward, economical, and high-throughput approach for observing thermal melts using commonly available real-time polymerase chain reaction instrumentation. However, quantitative analysis requires careful consideration of the dye-mediated reporting mechanism and the underlying thermodynamic model. We determine affinity constants by analysis of ligand-mediated shifts in melting-temperature midpoint values. Ligand affinity is determined in a ligand titration series from shifts in free energies of stability at a common reference temperature. Thermodynamic parameters are obtained by fitting the inverse first derivative of the experimental signal reporting on thermal denaturation with equations that incorporate linear or nonlinear baseline models. We apply these methods to fit protein melts monitored with SO that exhibit prominent nonlinear post-transition baselines. SO can perturb the equilibria on which it is reporting. We analyze cases in which the ligand binds to both the native and denatured state or to the native state only and cases in which protein:ligand stoichiometry needs to treated explicitly.

Thermodynamic analysis of a molecular chaperone binding to unfolded protein substrates.

Molecular chaperones are a highly diverse group of proteins that recognize and bind unfolded proteins to facilitate protein folding and prevent nonspecific protein aggregation. The mechanisms by which chaperones bind their protein substrates have been studied for decades. However, there are few reports about the affinity of molecular chaperones for their unfolded protein substrates. Thus, little is known about the relative binding affinities of different chaperones and about the relative binding affinities of chaperones for different unfolded protein substrates. Here we describe the application of SUPREX (stability of unpurified proteins from rates of H-D exchange), an H-D exchange and MALDI-based technique, in studying the binding interaction between the molecular chaperone Hsp33 and four different unfolded protein substrates, including citrate synthase, lactate dehydrogenase, malate dehydrogenase, and aldolase. The results of our studies suggest that the cooperativity of the Hsp33 folding-unfolding reaction increases upon binding with denatured protein substrates. This is consistent with the burial of significant hydrophobic surface area in Hsp33 when it interacts with its substrate proteins. The SUPREX-derived K(d) values for Hsp33 complexes with four different substrates were all found to be within the range of 3-300 nM.

Yeast screens identify the RNA polymerase II CTD and SPT5 as relevant targets of BRCA1 interaction.

BRCA1 has been implicated in numerous DNA repair pathways that maintain genome integrity, however the function responsible for its tumor suppressor activity in breast cancer remains obscure. To identify the most highly conserved of the many BRCA1 functions, we screened the evolutionarily distant eukaryote Saccharomyces cerevisiae for mutants that suppressed the G1 checkpoint arrest and lethality induced following heterologous BRCA1 expression. A genome-wide screen in the diploid deletion collection combined with a screen of ionizing radiation sensitive gene deletions identified mutants that permit growth in the presence of BRCA1. These genes delineate a metabolic mRNA pathway that temporally links transcription elongation (SPT4, SPT5, CTK1, DEF1) to nucleopore-mediated mRNA export (ASM4, MLP1, MLP2, NUP2, NUP53, NUP120, NUP133, NUP170, NUP188, POM34) and cytoplasmic mRNA decay at P-bodies (CCR4, DHH1). Strikingly, BRCA1 interacted with the phosphorylated RNA polymerase II (RNAPII) carboxy terminal domain (P-CTD), phosphorylated in the pattern specified by the CTDK-I kinase, to induce DEF1-dependent cleavage and accumulation of a RNAPII fragment containing the P-CTD. Significantly, breast cancer associated BRCT domain defects in BRCA1 that suppressed P-CTD cleavage and lethality in yeast also suppressed the physical interaction of BRCA1 with human SPT5 in breast epithelial cells, thus confirming SPT5 as a relevant target of BRCA1 interaction. Furthermore, enhanced P-CTD cleavage was observed in both yeast and human breast cells following UV-irradiation indicating a conserved eukaryotic damage response. Moreover, P-CTD cleavage in breast epithelial cells was BRCA1-dependent since damage-induced P-CTD cleavage was only observed in the mutant BRCA1 cell line HCC1937 following ectopic expression of wild type BRCA1. Finally, BRCA1, SPT5 and hyperphosphorylated RPB1 form a complex that was rapidly degraded following MMS treatment in wild type but not BRCA1 mutant breast cells. These results extend the mechanistic links between BRCA1 and transcriptional consequences in response to DNA damage and suggest an important role for RNAPII P-CTD cleavage in BRCA1-mediated cancer suppression.

Trafficking of G protein-coupled receptors.

G protein-coupled receptors (GPCRs) play an integral role in the signal transduction of an enormous array of biological phenomena, thereby serving to modulate at a molecular level almost all components of human biology. This role is nowhere more evident than in cardiovascular biology, where GPCRs regulate such core measures of cardiovascular function as heart rate, contractility, and vascular tone. GPCR/ligand interaction initiates signal transduction cascades, and requires the presence of the receptor at the plasma membrane. Plasma membrane localization is in turn a function of the delivery of a receptor to and removal from the cell surface, a concept defined most broadly as receptor trafficking. This review illuminates our current view of GPCR trafficking, particularly within the cardiovascular system, as well as highlights the recent and provocative finding that components of the GPCR trafficking machinery can facilitate GPCR signaling independent of G protein activation.

Smoothened signal transduction is promoted by G protein-coupled receptor kinase 2.

Deregulation of the Sonic hedgehog pathway has been implicated in an increasing number of human cancers. In this pathway, the seven-transmembrane (7TM) signaling protein Smoothened regulates cellular proliferation and differentiation through activation of the transcription factor Gli. The activity of mammalian Smoothened is controlled by three different hedgehog proteins, Indian, Desert, and Sonic hedgehog, through their interaction with the Smoothened inhibitor Patched. However, the mechanisms of signal transduction from Smoothened are poorly understood. We show that a kinase which regulates signaling by many "conventional" 7TM G-protein-coupled receptors, G protein-coupled receptor kinase 2 (GRK2), participates in Smoothened signaling. Expression of GRK2, but not catalytically inactive GRK2, synergizes with active Smoothened to mediate Gli-dependent transcription. Moreover, knockdown of endogenous GRK2 by short hairpin RNA (shRNA) significantly reduces signaling in response to the Smoothened agonist SAG and also inhibits signaling induced by an oncogenic Smoothened mutant, Smo M2. We find that GRK2 promotes the association between active Smoothened and beta-arrestin 2. Indeed, Gli-dependent signaling, mediated by coexpression of Smoothened and GRK2, is diminished by beta-arrestin 2 knockdown with shRNA. Together, these data suggest that GRK2 plays a positive role in Smoothened signaling, at least in part, through the promotion of an association between beta-arrestin 2 and Smoothened.

Purification, crystallization and preliminary X-ray diffraction studies of a complex between G protein-coupled receptor kinase 2 and Gbeta1gamma2.

G protein-coupled receptor kinase 2 (GRK2) phosphorylates activated G protein-coupled receptors (GPCRs), which ultimately leads to their desensitization and/or downregulation. The enzyme is recruited to the plasma membrane via the interaction of its carboxyl-terminal pleckstrin-homology (PH) domain with the beta and gamma subunits of heterotrimeric G proteins (Gbetagamma). An improved purification scheme for GRK2 has been developed, conditions under which GRK2 forms a complex with Gbeta(1)gamma(2) have been determined and the complex has been crystallized in CHAPS detergent micelles. Crystals of the GRK2-Gbetagamma complex belong to space group C2 and have unit-cell parameters a = 187.0, b = 72.1, c = 122.0 A, beta = 115.2 degrees. A complete data set has been collected to 3.2 A resolution with Cu Kalpha radiation.

Receptor and G betagamma isoform-specific interactions with G protein-coupled receptor kinases.

The G protein-coupled receptor (GPCR) kinases (GRKs) phosphorylate and desensitize agonist-occupied GPCRs. GRK2-mediated receptor phosphorylation is preceded by the agonist-dependent membrane association of this enzyme. Previous in vitro studies with purified proteins have suggested that this translocation may be mediated by the recruitment of GRK2 to the plasma membrane by its interaction with the free betagamma subunits of heterotrimeric G proteins (G betagamma). Here we demonstrate that this mechanism operates in intact cells and that specificity is imparted by the selective interaction of discrete pools of G betagamma with receptors and GRKs. Treatment of Cos-7 cells transiently overexpressing GRK2 with a beta-receptor agonist promotes a 3-fold increase in plasma membrane-associated GRK2. This translocation of GRK2 is inhibited by the carboxyl terminus of GRK2, a known G betagamma sequestrant. Furthermore, in cells overexpressing both GRK2 and G beta1 gamma2, activation of lysophosphatidic acid receptors leads to the rapid and transient formation of a GRK/G betagamma complex. That G betagamma specificity exists at the level of the GPCR and the GRK is indicated by the observation that a GRK2/G betagamma complex is formed after agonist occupancy of the lysophosphatidic acid and beta-adrenergic but not thrombin receptors. In contrast to GRK2, GRK3 forms a G betagamma complex after stimulation of all three GPCRs. This G betagamma binding specificity of the GRKs is also reflected at the level of the purified proteins. Thus the GRK2 carboxyl terminus binds G beta1 and G beta2 but not G beta3, while the GRK3 fusion protein binds all three G beta isoforms. This study provides a direct demonstration of a role for G betagamma in mediating the agonist-stimulated translocation of GRK2 and GRK3 in an intact cellular system and demonstrates isoform specificity in the interaction of these components.

Wee1-regulated apoptosis mediated by the crk adaptor protein in Xenopus egg extracts.

Many of the biochemical reactions of apoptotic cell death, including mitochondrial cytochrome c release and caspase activation, can be reconstituted in cell-free extracts derived from Xenopus eggs. In addition, because caspase activation does not occur until the egg extract has been incubated for several hours on the bench, upstream signaling processes occurring before full apoptosis are rendered accessible to biochemical manipulation. We reported previously that the adaptor protein Crk is required for apoptotic signaling in egg extracts (Evans, E.K., W. Lu, S.L. Strum, B.J. Mayer, and S. Kornbluth. 1997. EMBO (Eur. Mol. Biol. Organ.) J. 16:230-241). Moreover, we demonstrated that removal of Crk Src homology (SH)2 or SH3 interactors from the extracts prevented apoptosis. We now report the finding that the relevant Crk SH2-interacting protein, important for apoptotic signaling in the extract, is the well-known cell cycle regulator, Wee1. We have demonstrated a specific interaction between tyrosine-phosphorylated Wee1 and the Crk SH2 domain and have shown that recombinant Wee1 can restore apoptosis to an extract depleted of SH2 interactors. Moreover, exogenous Wee1 accelerated apoptosis in egg extracts, and this acceleration was largely dependent on the presence of endogenous Crk protein. As other Cdk inhibitors, such as roscovitine and Myt1, did not act like Wee1 to accelerate apoptosis, we propose that Wee1-Crk complexes signal in a novel apoptotic pathway, which may be unrelated to Wee1's role as a cell cycle regulator.

Suppression of conformational heterogeneity at a protein-protein interface.

Staphylococcal protein A (SpA) is an important virulence factor from Staphylococcus aureus responsible for the bacterium's evasion of the host immune system. SpA includes five small three-helix-bundle domains that can each bind with high affinity to many host proteins such as antibodies. The interaction between a SpA domain and the Fc fragment of IgG was partially elucidated previously in the crystal structure 1FC2. Although informative, the previous structure was not properly folded and left many substantial questions unanswered, such as a detailed description of the tertiary structure of SpA domains in complex with Fc and the structural changes that take place upon binding. Here we report the 2.3-Å structure of a fully folded SpA domain in complex with Fc. Our structure indicates that there are extensive structural rearrangements necessary for binding Fc, including a general reduction in SpA conformational heterogeneity, freezing out of polyrotameric interfacial residues, and displacement of a SpA side chain by an Fc side chain in a molecular-recognition pocket. Such a loss of conformational heterogeneity upon formation of the protein-protein interface may occur when SpA binds its multiple binding partners. Suppression of conformational heterogeneity may be an important structural paradigm in functionally plastic proteins.

A complex intronic enhancer regulates expression of the CFTR gene by direct interaction with the promoter.

Genes can maintain spatiotemporal expression patterns by long-range interactions between cis-acting elements. The cystic fibrosis transmembrane conductance regulator gene (CFTR) is expressed primarily in epithelial cells. An element located within a DNase I-hypersensitive site (DHS) 10 kb into the first intron was previously shown to augment CFTR promoter activity in a tissue-specific manner. Here, we reveal the mechanism by which this element influences CFTR transcription. We employed a high-resolution method of mapping DHS using tiled microarrays to accurately locate the intron 1 DHS. Transfection of promoter-reporter constructs demonstrated that the element displays classical tissue-specific enhancer properties and can independently recruit factors necessary for transcription initiation. In vitro DNase I footprinting analysis identified a protected region that corresponds to a conserved, predicted binding site for hepatocyte nuclear factor 1 (HNF1). We demonstrate by electromobility shift assays (EMSA) and chromatin immunoprecipitation (ChIP) that HNF1 binds to this element both in vitro and in vivo. Moreover, using chromosome conformation capture (3C) analysis, we show that this element interacts with the CFTR promoter in CFTR-expressing cells. These data provide the first insight into the three- dimensional (3D) structure of the CFTR locus and confirm the contribution of intronic cis-acting elements to the regulation of CFTR gene expression.

Computational Inference of Genome-Wide Protein-DNA Interactions Using High-Throughput Genomic Data

Transcriptional regulation has been studied intensively in recent decades. One important aspect of this regulation is the interaction between regulatory proteins, such as transcription factors (TF) and nucleosomes, and the genome. Different high-throughput techniques have been invented to map these interactions genome-wide, including ChIP-based methods (ChIP-chip, ChIP-seq, etc.), nuclease digestion methods (DNase-seq, MNase-seq, etc.), and others. However, a single experimental technique often only provides partial and noisy information about the whole picture of protein-DNA interactions. Therefore, the overarching goal of this dissertation is to provide computational developments for jointly modeling different experimental datasets to achieve a holistic inference on the protein-DNA interaction landscape.

We first present a computational framework that can incorporate the protein binding information in MNase-seq data into a thermodynamic model of protein-DNA interaction. We use a correlation-based objective function to model the MNase-seq data and a Markov chain Monte Carlo method to maximize the function. Our results show that the inferred protein-DNA interaction landscape is concordant with the MNase-seq data and provides a mechanistic explanation for the experimentally collected MNase-seq fragments. Our framework is flexible and can easily incorporate other data sources. To demonstrate this flexibility, we use prior distributions to integrate experimentally measured protein concentrations.

We also study the ability of DNase-seq data to position nucleosomes. Traditionally, DNase-seq has only been widely used to identify DNase hypersensitive sites, which tend to be open chromatin regulatory regions devoid of nucleosomes. We reveal for the first time that DNase-seq datasets also contain substantial information about nucleosome translational positioning, and that existing DNase-seq data can be used to infer nucleosome positions with high accuracy. We develop a Bayes-factor-based nucleosome scoring method to position nucleosomes using DNase-seq data. Our approach utilizes several effective strategies to extract nucleosome positioning signals from the noisy DNase-seq data, including jointly modeling data points across the nucleosome body and explicitly modeling the quadratic and oscillatory DNase I digestion pattern on nucleosomes. We show that our DNase-seq-based nucleosome map is highly consistent with previous high-resolution maps. We also show that the oscillatory DNase I digestion pattern is useful in revealing the nucleosome rotational context around TF binding sites.

Finally, we present a state-space model (SSM) for jointly modeling different kinds of genomic data to provide an accurate view of the protein-DNA interaction landscape. We also provide an efficient expectation-maximization algorithm to learn model parameters from data. We first show in simulation studies that the SSM can effectively recover underlying true protein binding configurations. We then apply the SSM to model real genomic data (both DNase-seq and MNase-seq data). Through incrementally increasing the types of genomic data in the SSM, we show that different data types can contribute complementary information for the inference of protein binding landscape and that the most accurate inference comes from modeling all available datasets.

This dissertation provides a foundation for future research by taking a step toward the genome-wide inference of protein-DNA interaction landscape through data integration.